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ABSTRACT 


A semi-inf ini te strip held rigidly on its short end is 
considered. Loads in the strip at infinity (far away from the 
fixed end) are prescribed. Integral transform technique is 
used to provide an exact formulation of the problem in terms 
of a singular integral equation. Stress singularity at the 
strip corner is obtained from the singular integral equation 
which is then solved numerically. Stresses along the rigid 
end are determined and the effect of the material properties 
on the stress intensity factor is presented. 

The method can also be applied to the problem of a lami- 
nate composite with a flat inclusion normal to the interfaces. 
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INTRODUCTION 


In the past there has been considerable interest in the 
strip problems from both mechanics and mathematical points of 
view. A few good solutions exist [1-6]; however, none of these 
methods appears to provide a solution which can directly give 
the correct behavior of the stresses at the corners and does 
not present any convergence difficulties in the evaluation of 
the contact stresses. In [1-4], the stress singularities at 
the corners were ignored. Benthem [5] used the Laplace trans- 
form technique to formulate the problem but required a knowl- 
edge of stress singularity by alternate means [7], In [6], 
Vorovich and Kopasenko considered a problem of a semi - i nf i ni te 
strip loaded symmetrically on its longitudinal sides and re- 
duced the problem to a singular integral equation. Again 
stress singularities were obtained by separate analysis and 
the integral equation was solved using a very cumbersome method. 
With the recent advances in the analysis and solution of such 
singular integral equations [8], it seems possible to extract 
the corner si ngul ari ties and solve the derived integral equa- 
tion in [6], The method presented here reduces the problem to 
a singular integral equation which contains the stress singu- 
larity and can be solved numerically giving very satisfactory 
resul ts . 

To solve the problem of a strip with fixed end we will 
first consider the plane elastostatic problem for a flat inclu- 
sion situated centrally in an infinite strip with sides free 
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of tractions (Figure 1). A similar problem with the inclusion 
lying parallel to the interface was considered in [9] where 
stress field around the inclusion tip was developed. Since 
the problem considered here is symmetric with respect to the 
axis containing the inclusion, a rigid inclusion has the same 
effect as an inextensible inclusion with zero bending rigidity. 
The case of an elastic inclusion can also be treated without 
any added difficulty [9]. In this paper, we will concentrate 
on the case when the inclusion extends to the surfaces, thus 
reducing the problem to that of a semi - i nf i ni te strip with 
fixed end. Numerical results will be presented only for the 
strip problem. 

The method can easily be applied to the case where the 
strip containing the inclusion is bonded to two similar half 
planes with different elastic properties than the strip. This 
problem makes the equations and the algebra much more involved, 
but presents no conceptual difficulty. 

FORMULATION OF PROBLEM 

Consider a strip of width 2h in plane strain condition 
with a rigid or inextensible inclusion of length 2a (a s h) 
located centrally along the x-axis (Figure la). The shear 
modulus and Poisson's ratio of the strip are y and v respec- 
tively. Let the strip be uniformly loaded at both ends far 
away from the inclusion. Figure lb shows the main problem of 
interest which is obtained from Figure la by extending the 
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inclusion to the surfaces x = ±h, or letting a = h. Note that 
the chosen problem has symmetry with respect to both x and y 
axes . 

The problem described above can be recovered by the super- 
position of two subproblems I and II as shown in Figure 2. In 
I, we simply have a uniformly loaded strip without any inclu- 
sion. Solution of I is given as 

“yy’tx.y) * °yy 1 ( X » “ ) = P Q 

c-xy^x.y) = o xx '(x,y) = 0 

( 1 ) 

u,(x,0) = -e 0 x ; p 0 

v j ( x , 0 ) = 0 

where k = 3-4v for plane strain and k = (3-v)/(l+v) for plane 
stress. 

Problem II is the disturbance problem where the input 
function must be the displacement in x-direction at y = 0 plane, 
equal to the negative of that in I. Hence the boundary condi- 
tions for II become 

°xx <h ’ y) = a xy 

v ( x ,0) = 0 , 

u(x,0) = e Q x , 

a xy ( x ’0) = 0 , 


(h,y) = 0 , 
| x | < h , 


( 2 ) 


I x | < a , 
a < | x | < h , 
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Due to the symmetry of the problem, it is sufficient to 
consider only one quarter of the problem. Note that y = 0, | x | < a 
is a singular surface across which the displacements are con- 
tinuous and the stress vector suffers a discontinuity. The 
displacement and stress fields for the strip can be written as 
a superposition of two transform solutions. One is the solu- 
tion for an infinite strip (|x)<h, |y|<°°) in plane strain with 
symmetry about x=0 and y=0 planes, and the other is the half 
plane solution with x=0 as the plane of symmetry. Expressing 
the solution as 

CO _ _ 

u(x,y) = - f / {kf(n) - g (n ) ] si nh(nx) 

+ xg(n)cosh(rix) Icospy dn ~ ~ J (<-£y)e"^ y sin£xd£ 

o ^ 

v(x,y) = | / {l[f(n) + g(n)]cosh(nx) 

2 °° r v 

+ xg(n)sinh(nx)}sinny dn + - / yd> ( C ) e ^cos^x d E, 


a xx (x ’ y) 


2y 


- - / [f (n)cosh(nx) + nxg(n)sinh(nx)]cosny dn 
o 


- | 7 - ?y]e' 5y cosCx dK 

o 


a (x,y) 9 oo 

- 11 = | / {[f(n) + 2g(n)]cosh(nx) 


(4) 


2y 


+ nxg(n)sinh(nx)}cosnydn + / 4>(£)[^r" " ^ y cos£xd£ 
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a (x,y) 2 


/ {[f (n) + g(n)]sinh(nx) 


+ nxg(n)cosh(rix)}sinriydri + ~ f 4>U)C^y^ - £y]e~^ y si n£xd£ 


it may be seen that this solution identically satisfies the 
third condition v(x,0) = 0 of (2). The three unknowns f, g and 
4) must then be determined by using the first two conditions of 
(2) and the mixed boundary conditions (3). The first two con- 
ditions of (2) may be written as 

f(n)cosh(nh) + nhg(n)sinh(nh) 


| 7 Sit# [EJl ♦ -^cosChdC 
0 n +€ n +5 


f (n)sinh(nh) + g ( n ) [sinh(nh) + oh cosh(nh) 


o n +£ n +£ 


Mixed boundary conditions (3) can be expressed as follows: 


iu(x ,i 

9x 


- - / (K-Cy)4>U)e ^ y cos£xd£ 

o 

1 i m 
y+0 + 


O 00 q 

- £ / { [f ( n ) - g(n)]cosh(nx) + oxg(o)sinh(nx)}dn 


x < a 


5 yu ( x »0 + ) ^ + i °? , , . , 

-^5 = ^r~ ! <K?)sin£xd£ = 0 , a < f x | < 
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Note that displacement derivative is used in (6) instead of 
displacement in order to have a dimensional consistency in (6) 
and ( 7 ) . 


A new unknown function G(x) is now defined as the shear 
stress along the line y = +0 which is to be distinguished from 
the shear stress at y = -0 as the stress vector suffers a dis- 
continuity across this boundary. Due to symmetry, the two 
values should be numerically equal and negative of each other. 
The function G(x) is written as 

CT xy(x ,0+ ) = G ( x ) » M < h (8) 

From (7), G(x) =0, a < Jx| < h. Inverting the integral ob- 
tained from (7) and (8) we have 

l ® 

4>U) = - I G(t)sin£tdt (9) 

Solving (5) simultaneously we obtain 


i D, (n) [si nh(nh) + nh cosh(nh)] - D 0 (n)nh sinh(nh) 

i f(n) ■ -J * 

sinh(2nh) + 2qh 

( 10 ) 

-D, (n)sinh(nh) + D 9 (n)cosh(nh) 

7 g(n) = — 1 ? 

c Sinh(2nh) + 2nh 


where 


D,(n) • - f 1 ^4 [SJ1 ♦ ]cos 5 hdc 

1 71 0 n +5 “ 


2. e 2 
n +5 


(ID 


° 2 {n) = - | 7 D |^- ]sin?hde 

c o n +5 n H 
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In order to reduce (6) to a singular integral equation 
f and g must be substituted in it from (10) and then (9) should 
be used to obtain the equation in terms of G(x). In this sym- 
metric problem shear stress is an odd function in x; therefore, 
G(x) = -G(-x). Using this property, the first integral in (6) 
when combined with (9) yields a Cauchy kernel: 


lim / G(t)dt / e"^ y s i n£{ t-x) d£ = / k( . p . dt (12) 

y+0 + -a o -a t_x 

Now substituting <j> from (9) into (11) and changing the order 
of integration gives 

2u(k+ 1 )D 1 (n) 


= | f G(t)dt / ]s 1 n€(h-t)d 5 

* -a o n +S d n+5 

( 13 ) 

2y ( k+ 1 ) D 2 ( n ) 

= - | / G(t)dt / [^- - -P-T ]cosC(h-t)dC 

-a o n +? n +5 

From the tables of Fourier Transform in [10] we have 


cos ?z 

2 , r 2 
n +£ 


d? = 


TT 

2n 


-nz 


l ~tn sinczd5 = ze ' nz 

0 (n+S 2 ) 


( 14 ) 


By differentiating (14) with respect to z we obtain 
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7 "o— y sinCzdC = £ e" nz 
o r/+r 2 

(15) 

00 r 2 

/ cos^zd? = J- [1 - nz]e' nz 

0 (nV> 

Using (14) and (15) and replacing z by h-t, (13) can now be 
written as 

2u(k+1 )D 1 (n) = / G(t)[^- + n(h-t)]e“ n ^ h " t ^dt 

(16) 

2)a ( K+l ) D 2 (n ) > - I G(t)[E|i + T ) (h-t)]e' n(h ‘ t) dt 

The singular integral equation obtained from (6) using relations 
(12), (10) and (16) can’ now be expressed as 

k / dt + / G(t)K(t,x)dt = - y ( k+ 1 ) ire , 

-a -a 

I x | < a (17) 

where 

K(t,x) = / k(t,x,n)e“ T1 ^ h ' t ^dri 
o 

e -nh r 

k(t,x,n) = ottu 2T7F cosh (nx) {1 

.1 + 4nhe" 2nn - e 4nn L 

+ [k + 2ri(h-t)][K - 2 + 2nh] + 2e 2ri ^(l - k + qt)} 

- 2qxsi nh(qx ) {k + 2q(h-t) - e” 2nll }J (18) 

The kernel K(t,x) is bounded for all values of t and x in 
(-a, a) if a < h. The integral equation (17) should be solved 


- 9 - 



subject to the following equilibrium condition: 
a 

/ G ( t ) dt = 0 (19) 

-a 

For a < h, the solution of (17) is rather straightforward [8], 
and will not be considered in this paper. The actual problem 
of interest (Figure lb) is the case when a = h. For this case 
the kernel K(t,x) is no longer bounded for all values of t and 
x, and contains point singularities at t = h and x = ±h. These 
singularities can be extracted by using the asymptotic value 
of the integrand k(t,x,n) as n + °» in (18). Let 

K (t,x) = / k (t,x,n)e“ n ^ h " t )dn (20) 

s o 

where K s (t,x) is the singular part of the kernel. From (18) 
it follows that 

2 

Mt.x.n) = 2e" 7lh ^cosh(nx){ ^ lc - ^ 1 + nUh + (k-2) (h-t)] 

+ 2n 2 h(h-tj } - nx sinh(nx){< + 2n(h-t)} (21) 

Using the following result [10] 

J ri m e -n(2h-t) C sinh(nx)'! dr) = d^_ J e -n(2h-t) ( si nh(nx)'| dn 
o v,cosh(nx)J dt m o Lcosh(Tix)J 

x 

2h-t 

the singular kernel K (t,x) now becomes 


= r i 

dt m (2h-t) 2 - x 2 
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K s (t,x) 


(k- 1 ) 2 (2h-t) 
[(2h-t) 2 - x 2 ] 


+ 2 + ( K ' 2 )( h - t )>{(2h-t) 2 + x 2 } 

[(2h-t) 2 - x 2 ] L 


- 2 K (2h-t)x 2 + 40wt} — {h(2h-t)[(2h-t) 2 + 3x 2 ] 

( 2h - 1 / - x^ 

- x 2 [3(2h-t) 2 + x 2 ]] (23) 

Also, we observe that k(t,x,n) + 00 near n = 0. This can be 
easily handled by isolating the value of the integral near 
n = 0. We have 


K (t,x) = / k (t,x ,n)e‘ n ^ h_t ^dn 
o 

where 


k 0 (t,x,n) 


(3-k) ( 1 -k) 
8nh 


Therefore , 

* 0 (t.x) = - ln(h-t) 


(24) 


(25) 


Note that K (t,x) is independent of x and has a weak singular- 
ity in the sense that it is square integrable, hence does not 
need any special attention. 


To analyze the behavior of the unknown function G ( t ) near 
the end points, dominant part of the equation consisting of 
the Cauchy kernel and the singular kernel K g (t,x) must be con- 
sidered. We can express these dominant terms as 


- 11 - 



? j te + 7 <* 2 - 3 ♦ 12("-x)Sx - 4(h - x)2 ^7> 7F^ 

+ 7 {k2 ' 3 - 12 < h+x >33T - 4 < h+x > 2 7T> 7TfW] dt ' 

n Y — J 


= - y ( k+ 1 )e 0 + A(x) , 


x | < h 


(26) 


where A(x) is a bounded function containing the terms coming 
from the Fredholm kernel in (17), i.e., 
h 

A ( x ) = - / [ K( t ,x) - K ( t ,x) ]G( t )d t 
-h s 

G(t) is assumed to have integrable singularities at t = ±h and, 
following [11], may be expressed as 


ftf<h 

(h 2 -t 2 ) (t-h) a (t+h) a 


(27) 


where 0 < Re(a) < 1 and H ( t ) satisfies a Holder condition in 
the closed interval -h < t < h. Procedure for determining a 
requires studying (26) as in [11, Chapter 4] and has been pre- 
sented in detail in [8]. 


Considering the following sectionally holomorphic function 


*<z) - I /f^dt . 1 / 


H( t)e 


-h (t-h) a (t+h) a (t-z) 


(28) 


according to [8, Chapter 4] 


<f>(z) 


H(-h) e^ 


,a sirnra 


(z+h)< 


^ 5 a + 

( 2h) a sinTra(z-h) a 0 


(29) 
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where <J> Q (z) is bounded everywhere-' except possibly at the end 
points ±h where it may have the following behavior: 

H (±h) 

1 4> 0 (z) I < — — , Re ( a Q ) < Re ( a) (30) 

(z±h) 0 

Now taking the limit of (29) as z -> x , 2h+x and 2h-x, we find 


4>(x ) = 

H ( -h ) cot-rra 
(2h) a (h+x) a 

♦ ♦*(*) . 1 x 1 < h 

(2h) a (h-x) a 

<f>(2h+x) 

H ( h ) 

1 

♦ /(x) » h < 2h+x < 3h 

(2h) a sin7ra (h+x) a 

<p( 2h-x) 

H ( h ) 

1 + 

4> ^ * ( x ) * h < 2h-x < 3h 


... /U „\a 


(2h)°siniTa (h-x) 


(31) 

it ★ "it 

where around the end points <p , <p - j , and <p ^ have behavior 
similar to that of Q ( z ) . 

Substituting (31) into (26) we obtain 


[kH (- h )cosira + ^ - ^ -{(k^- 3) + 12a - 4a(a+l)}] 

(2h) a sin-rra(h+x) a 2 


H(h) 

(2h) a sinTra(h-x) a 


[ kc o s na - ^-{(k^- 3) + 12a - 4a(a+l)}] 


where P(x) contains all the bounded functions. 

Since G(t) is an odd function of t, hence H(t) = - H(-t), 
and since H(±h) f 0, equation (32) can only be satisfied if 

2KCOS-rra = + 1 - 4(a-l)^ (33) 
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which is the characteristic equation to determine a. The equa- 
tion is identical to the one derived in [7] for a wedge with 
stress-displacement boundary conditions. It depends only on 
the Poisson's ratio of the strip and gives a real value of a 
for any material, i.e., for 0 < v < 0.5. If a < h as in Figure 
la, only singular kernel is a Cauchy kernel and the character- 
istic equation becomes 

cotiTa = 0 , a = 1/2 (34) 

which is the well known singularity at the tip of a flat inclu- 
sion [9]. For the strip problem, the singular integral equa- 
tion can be written as 

h 

dt + ^ G(t)[K $ (t,x) + K Q (t,x)]dt 

+ ^ G( t) Kp( t ,x)dt = - y(K + l)Tre 0 , |x| < h (35) 

Kp( t ,x) = 7 [k(t,x,n) - k 0o (t,x,n) - k o (t,x,n)]e" r| ^ h " t ^dt 

K s (t,x), K Q (t,x), k co ( t ,x ,n ) and k Q (t,x,n) are given by equa- 
tions (23), (25), (21) and (24) respectively. Kp(t,x) is a 
Fredholm kernel for a < h. 

SOLUTION OF THE INTEGRAL EQUATION 

To solve (35) we first normalize the dimensions with 
respect to h by the following transformations: 


ii 

K i 


G(t) 

t-x 
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T = £ > y = £ , G(t) = G(hx) = <J>(t) (36) 

Hence (35) can be expressed as 

1 i h 

| 4>(r) [^Ty + £ '{K s (hT,hy) + K 0 (hx,hy) + K p ( hx , hy ) }]dx 


_ _ y(K+1 )tt£ 0 
K 

and (27) becomes 


I y I < 1 


(37) 


4>( x) 


d-x 2 ) a 


(38) 


where a is given by (33). Equation (37) can now be solved by 
using Gauss-Jacobi Integration formula. The method has been 
previously used and described in [8], Equation (37) must be 
solved subject to the additional condition (19). We obtain a 
set of NxN simultaneous algebraic equations given as 


N 

l 

j = l 


V (T j )[ 


T j' y i 


h 

K 


{ K $ ( hxj ,hy i ) + K Q ( hxj ,hy i ) 


+ K p (hx j ,h yi )}] 


p(k+1 ) tte q 

K 


N 

l A • ip ( x - ) = 0 
j = I J J 


(39) 


where, noting that (l-x^) a is the weight function of Jacobi 
polynomials " a,_a ^ (x) , the constants x j and y^ were shown 
to be the roots of the following equations [9]: 
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P N { ' a, ' a)(x j) = 0 ’ (J = 

P N _l ( 1 ' a, 1 “ a) (yi) = 0 , (i = 1 N-l) 

and A.'s are the corresponding weighting constants [9]. iMt.) 

J J 

are numerically computed from (39). The shear stress a (x,0) 

xy 

can then be expressed as 

h 2 a ^(r-) 

a xv (x,0) = G(x) = V— , | x | < h (40) 

* J 9 9 a 

(h 2 -x 2 ) 

NORMAL STRESS a yy (x,0) AND STRESS INTENSITY FACTOR 

After solving for the shear stress for the disturbance 
problem, the only remaining important quantity of interest is 
the normal stress a yy (x,0). To evaluate a yy , we start with 
the fourth equation in (4) and use (9), (10), (12) and (16). 
For y = 0 we find 

ff yy (x,0) = §T 1 dt + £TT l _[ G < t )t K io (t ’ x) 

+ K 1$ (t,x) + K-j F ( t,x)]dt , | x J < h (41) 

where 

K-jp(t,x) = 

7 C k i ( t ,x ,q) - k loo (t,x,n) - k 10 (t,x,n)]e" n ^ h " t ^dn 
o 

and 
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K 1Q (t,x) = " log(h-t) 


K 1 s (t,x) 


( 1 -3k) ( 2h- 1 ) 
(2h-t) 2 - x 2 


+ 1 

[ ( 2h- t) 2 - 

- 2i<(2h-t)x 2 

- x 2 [3( 2h-t) 2 


' 2 ] 


2 { Kh - 3(h-t)H(2h-t) 2 + x 2 } 


+ {h(2h-t) [(2h-t) 2 + 3x 2 ] 

(2h-tp - x ^ 


+ x 2 ]} 


(42) 


,-nh 


kl(t ’ X,n) = 1 + 4nhe- Z ^ - e" 4 ^ L' 


cosh(nx) {1 


[ tc + 2ri(h-t)][3 - 2rih] + e 2ri ^(3 - k + 2nt)} 


- 2nxsinh (nx) {k + 2ri(h-t) - e 2ti ^}J 

kioo(t» x »n) = e _nh £cosh(nx){l - [k + 2n(h-t)](3 - 2n h ) 

- 2nxsinh(nx){K + 2ri(h-t)}J 

k 10 (t,x,n) = ^ 

To determine the behavior of a yy( x >0) near the corner points, 
we must consider the dominant part of the equation (41) which 
can be written as 

( K+ 1)o yy (x,0) = 1 J 6(t)[^l 

2 

+ ‘{-(3k+5) + 2 ( k+7 ) (h-x)^ - 4( h-x) 2 } 2h ] - t — 

d x 

2 

+ {-(3k+ 5) - 2(<+7)(h+x)^- - 4( h+x) 2 } 2h- ' t+x ^ dt 

d x 

(43) 
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Substituting (27) in (43) and using the relations (31), the 
dominant part of the normal stress becomes 


( K+ 1)a (x,0) = 1 [(k-1 ) (cosirot + 1) - 2 ( k+1 ) (a-1 ) 

yy ( 2h ) a s i nira 

+ 4(a-l) 2 ][ ] , x - ±h (44) 

(h+x) a (h-x) a 

Stress intensity factors K-| , K 2 can be defined as 

K ] = lim /2 (h-x) a cr (x ,0) 
x-*h yy 

(45) 

K 2 = lim ( h-x) a a (x ,0) 
x+h xy 


Using (27) and (44), these can be rewritten as 

K = i^2_HihJ [{<“1 ) (cos Trot + 1) - 2 ( k+ 1 ) (a-1 ) 

(k+ 1 ) (2h) a simra 


+ 4(a-ir] 


K = n Mil)_ 
" (2h) a 


(46) 


NUMERICAL RESULTS AND DISCUSSION 

The total solution of the problem shown in Figure lb is 
now the sum of the two problems I and II. Hence, 


and 


a yy T(x ’ 0) = a yy I(x,0) + a yy (x,0) = p o + a yy (x>0) 


a xy T (x ,0) = a xy I ( x>0 ^ + a xy^ x, °^ = a xy^ x, °^ 


(47) 
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Since the problem of main interest here is the semi -i nf i ni te 
strip constrained at y = 0, the results only for the case of 
(a = h) are presented. Figures 3 and 4 show the variations 
of normal and shear stresses, respectively, along the fixed 
end for various values of the Poisson's ratio. When the Pois- 
son's ratio is zero, we do not have a disturbance problem and 
the solution of the total problem is identical to that of 
problem I (Figure 2). Thus, for v = 0 


o xy T (x,0) = 0 
° y y T ( x >°) ' P 0 


(48) 


As the Poisson's ratio of the strip increases (to a maximum 
value 0.5), effect of the disturbance problem (II in Figure 2) 
increases as shown in these figures. Also, for higher value 
of the Poisson's ratio, a larger value of the power a of the 
stress singularity is obtained. This effect is depicted clear- 
ly from graphs 3 and 4 by the behavior of stresses near x = h. 


Figure 5 shows the variation of the stress intensity 
factor (as defined in (46)) with respect to the Poisson's 
ratio of the strip. As seen from (46), the stress intensity 
factor K-j depends on «£ and their ratio l^/Ki is a function 
of the power of the stress singularity a and the Poisson's 
ratio. From (46) 

_ 2 _ ( k+ 1 ) s i mra ^ 

^1 [ ( k- 1 ) ( costo + 1) - 2(k+1 ) (ot-1 ) + 4(a-l)^] 
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Figure 5 also shows a variation of this ratio K 2 /K.j with re- 
spect to v. Negative sign in (49) comes from the fact that 
we have defined the stress intensity factors near x = h (45) 
where we have negative shear stress and positive normal stress. 
The physical significance of the ratio K 2 /K-| is clear in the 
case of an elastic strip pressing against a much stiffer body 
(i.e., p Q is negative). In this problem if the coefficient 
of friction, f, is greater than K 2 /K-|, it may be assumed that 
no sliding would occur between the strip and the adjoining 
body, i.e., the contact condition is that of perfect adhesion, 
and the solution given in this paper would be valid. If 
f < K^/K-j , the problem becomes that of an elastic punch on a 
rigid half space with friction. From Figure 5 it is seen that, 
in testing whether the end condition is that of perfect adhe- 
sion or sliding for a given strip under compression, as a first 
approximation one may assume that K 2 /K-j = v. 
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Figure 3. Shear Stress vs. Poisson's Ratio for 
the Semi - Inf ini te Strip. 



Figure 4. Normal Stress vs. Poisson’s Ratio for 
the Semi - Inf i ni te Strip. 






